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Аннотация 

Введение. Основой для исследования, анализа и математической оптимизации любого химического процесса 
является адекватная математическая модель, учитывающая кинетику объекта. Кинетический анализ 
в химической технологии является важной задачей, поскольку позволяет оптимизировать процессы синтеза 
и прогнозировать их эффективность. Многие химические процессы включают в себя несколько стадийных 
реакций. Для успешного проектирования и оптимизации необходима математическая модель, которая 
описывает каждую стадию. Создание такой модели вручную может быть трудоемким и затратным процессом, 
требующим обработки большого объема информации. Современный уровень автоматизации позволяет 
ускорить получение математического описания кинетики многостадийных реакций. В этом случае значительно 
упрощается работа с данными и уменьшается вероятность совершения ошибок. Полученная математическая 
модель может быть применена для последующего анализа и оптимизации процесса. В работе рассмотрена 
промышленная реакция каталитического риформинга бензина, занимающая важное место в современной схеме 
переработки нефти, поскольку является источником высокооктановых компонентов товарных бензинов 
и индивидуальных ароматических углеводородов. Данный процесс характеризуется участием в нем большого 
числа (до 300) различных углеводородов, изменением числа молей и неизотермичностью. Математическое 
моделирование таких процессов предполагает детализацию стадий до необходимого уровня. Рассмотрена 
детализация до 173 стадий. В такой постановке задачи автоматизация формирования математического описания 
кинетики для каталитического риформинга бензинов ранее не проводилась. Поэтому целью представленной 
работы явилась реализация эффективных численных методов и алгоритмов для автоматизации формирования 
математической модели с учётом кинетики, термодинамики и изменения числа молей. 

Материалы и методы. Математическое описание кинетики многостадийных реакций разрабатывается 
на основе закона действующих масс. Значения кинетических параметров взяты из литературных источников. 
Решение прямой задачи кинетики проводилось с применением следующих алгоритмов: метод Гира, Рунге- 
Кутты 4 порядка и метод зс1ру.одеш языка Руфоп. Концепция автоматизации реализована с помощью 
методологии ШЕЕЙ. Программное обеспечение написано на языке программирования Ру®оп. 

Результаты исследования. Создано новое программное обеспечение для автоматизации процесса 
формирования математической модели с учетом кинетики, термодинамики и учета объема реакционной смеси. 
Приведены результаты работы программы на примере каталитического риформинга бензина. Реализована 
возможность учета в модели промежуточного подогрева смеси в каскаде реакторов. Получены численные 
значения изменения температуры, соответствующие промышленным данным. 

Обсуждение и заключение. Результаты, полученные при моделировании химических превращений в каскаде 
реакторов каталитического риформинга бензина, подтвердили экзотермический характер реакции. 
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Разработанный программный продукт позволяет отобразить изменения концентраций веществ реакции, а также 
изменение температуры в реакторе и может быть использован в научных исследованиях организаций ДлЯ 
анализа многостадийных каталитических процессов. Результаты моделирования кинетики реакции будут 
использоваться в последующей оптимизации условий проведения процесса на производстве. 
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Введение. Изучение кинетики многоступенчатых химических процессов имеет теоретическое 
и практическое значение, поскольку позволяет раскрыть механизм реакций и осуществить контроль. При этом 
становится возможным ускорять целевые или замедлять побочные химические реакции. Кинетический анализ, 
будучи важной задачей химической технологии, способствует оптимизации процессов синтеза 
и прогнозированию их эффективности [1-3]. Многие химические процессы являются многостадийными 
реакциями. Для их эффективного проектирования и оптимизации необходимо иметь математическую модель, 
описывающую все стадии [4,5]. Ручное формирование такой модели может быть очень трудоемким 
и затратным процессом, так как подразумевает под собой обработку большого объема информации [6]. 

Объектом исследования является промышленный процесс каталитического риформинга бензина. Целевой 
продукт реакции — риформат — позволяет получать основную массу товарного бензина. Остальные 
составляющие являются продуктами изомеризации, каталитического крекинга и алкилирования. Большинство 
нефтеперерабатывающих заводов в России и мире используют установки каталитического риформинга [7, 8]. 
Каталитический риформинг бензина позволяет получать ароматические углеводороды (ксилол, толуол, бензол), 
которые затем могут быть использованы в нефтехимической отрасли [9, 10]. Однако для товарного бензина 
присутствие в продукте ароматических углеводородов ограничивается по экологическим критериям. Поэтому 
детализированная сложная кинетика, показывающая постадийное образование компонент, позволит управлять 
условиями процесса для контроля выхода продукта и ограничения содержания ароматических углеводородов. 
Процесс каталитического риформинга бензина проводится в реакторном блоке, состоящем из трех 
последовательных реакторов, между которыми установлены печи для нагрева смеси. Это является 
необходимостью, так как данная реакция — эндотермическая. При математическом описании реакторного 
блока необходимо, наряду с изменениями концентраций компонент от условного времени контакта 
катализатора и смеси реакции, учитывать изменение температуры и объема реакционной смеси. 

Автоматизация моделирования таких сложных процессов позволяет значительно ускорить получение 
математического описания кинетики многостадийных реакций, снизить вероятность ошибок и упростить 
процесс работы с данными [11]. В целом автоматизация формирования математического описания для 
многостадийных химических реакций является важным инструментом для эффективного управления 
процессами в различных отраслях промышленности, а также для достижения более точной и эффективной 
реализации химических реакций в лаборатории. 

Таким образом, целью исследования являлась разработка программного продукта, помогающего 
автоматизировать формирование математического описания многостадийных химических реакций с учётом 
кинетики, термодинамики и изменения объёма реакционной смеси, а также получение численного решения 
составленной математической модели. 

Для достижения цели исследования в данной работе решались следующие задачи: 

— разработка программного модуля формирования математического описания кинетики с использованием 
языка программирования Руфоп; 

— разработка программного модуля численного решения прямой задачи согласно составленному 
математическому описанию с использованием языка программирования Ру®оп; 

— получение результатов работы программного продукта для реакции каталитического риформинга бензина. 

Материалы и методы. Математически указанный процесс может быть описан системой обыкновенных 
нелинейных дифференциальных уравнений с начальными данными, которая и является математической 
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моделью. Учитывая, что каталитический риформинг бензина является неизотермической реакцией 
с изменением объема реакционной смеси, модель многоступенчатой химической кинетики может быть 
рассмотрена в следующем контексте: 
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уй 
н,(Т)=АН, (298) [С„(Т)ат; (6) 
298 
С„(Г)=а, +В Т+сТ +а. 1"; (7) 


с начальными условиями при т= 0, у, (0) = у’, Т = 763К, О =1.014; те [о,г | ПИ. 


В данной формуле реакционной кинетики имеется ряд переменных и коэффициентов, которые 


определяют концентрацию реагентов, время контакта и скорость реакции. Конкретно О иу, (концентрации 


реагентов реакции) могут быть выражены в моль/л или в долях, а переменная т, показывающая условное 
время контакта реакционной смеси с катализатором, в кгхмин/моль или кг катализатора. В формуле также 
задействованы показатели: количество стадий М (173 стадии) и количество веществ № (38 веществ), а 


Также стехиометрическая матрица У; . (. 


в = это все отрицательные элементы матрицы У; на пересечении 


7-ой строки и /-го столбца, а В, — все положительные элементы матрицы У, на пересечении {-ой строки и 
7-го столбца. Скорость /-ой стадии реакции у, может быть выражена в 1/мин или в моль/(кгхмин), а 
константы скоростей стадий А.К, — в 1/мин. Кроме этого, в формуле указаны предэкспоненциальные 
множители в. м выраженные в 1/мин, энергии активации прямой и обратной реакции Е,, Е), 
выраженные в ккал/моль, газовая постоянная К и температура Т, измеряемая в Кельвинах. 


Продолжительность реакции 1* определяется временем, которое требуется для её завершения в минутах, а 


энтальпия образования 1-го компонента АН, (Т) измеряется в Дж/моль. Удельную теплоемкость С’, (Т) 7-го 


компонента, измеряемую в Дж/(мольхК), можно выразить через коэффициенты температурной 
зависимости теплоемкости данного компонента, которые обозначены буквами а,Б,с,,а, [12]. 1 — 


порядковый номер вещества, участвующего в реакции. 

При моделировании химических превращений каталитического риформинга бензина необходимо 
учитывать болышое число индивидуальных углеводородов. Поэтому в работе применяются групповые 
компоненты в модели по выделению типов углеводородов и числу атомов углерода в структуре молекулы: 
А; — ароматические углеводороды; АСР, — пятичленные нафтены; АСН; — шестичленные нафтены; 1Р, — 
изомеры парафинов; пР; — нормальные парафины. Здесь / — количество атомов углерода в структуре 
молекулы. В таблице | приведены усредненные значения коэффициентов температурной зависимости 
теплоемкости групповых компонент, которые будут применяться в расчетах. 
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Таблица 1 
Коэффициенты температурной зависимости групповых компонент 
| ре а, Ь, с; а, 
углеводородов 
1 ПР! 34,942 —0,03996 0,19184 —0,153 
2 пПР› 28,146 0,043447 0,18946 —0,1908 
3 ПР: 28,277 0,116 0,19597 —0,2327 
4 пР4 20,56 0,2815 —0,013143 —0,09457 
5 ТР. 6,772 0,34147 —0,10271 —0,03685 
6 ПР; 26,671 0,32324 0,04282 —0,1664 
7 1Р —0,881 0,47498 —0,24797 0,06751 
8 пРе 25,924 0,41927 —0,012491 —0,1592 
9 ТР —7,123 0,58327 —0,30338 0,06802 
10 ПР7 26,984 0,50387 —0,04748 —0,1684 
И 1Рт 19,245 0,55072 —0,14055 —0,08248 
12 ПР 29,053 0,58016 —0,057103 —0,1955 
13 1Рз —3,367 0,75824 —0,38216 0,05736 
14 пПРо 29,687 0,66821 —0,096492 —0,2001 
15 1Рэ 68,581 0,44754 0,31908 —0,5118 
16 пРо 31,78 0,74489 —0,10945 —0,2267 
17 Ро —46,17 1,108 —0,70316 0,1787 
18 пПРи 125,21 0,31401 0,79137 —0,9141 
19 Ри —8,791 1,0548 —0,5778 0,1192 
20 АСНб 13,783 0,20742 0,53682 —0,6301 
21 АСНУ? 4,296 0,42716 0,21058 —0,3999 
22 АСНв —51,866 0,78827 —0,35255 —0,006855 
23 АСНо —120,89 1,2728 —1,0794 0,4035 
24 АСНуо 90,421 0,23264 0,94595 —1,057 
25 АСНи —65,48 1,1809 —0,68379 0,1458 
26 АСР —9,939 0,42528 0,012521 —0,1886 
27 АСР7 —28,514 0,58607 —0,094379 —0,1644 
28 АСР& —69,713 0.92602 —0,62526 0,2049 
29 АСРо —116,73 1,3097 —1,2439 0,5292 
30 АСР —67,341 1,0922 —0,704 0,1906 
31 АСРи —68,23 1,187 —0,7575 0,2018 
32 Аб —31,368 0,4746 —0,31137 0,08524 
33 А? —24,097 0,52187 —0,29827 0,06122 
34 Аз —17,36 0,5647 —0,26293 0,01122 
35 Ао —10,933 0,64349 —0,27829 —0,01443 
36 Аю —24,187 0,79716 —0,48265 0,1341 
37 Ап —26,717 0,91044 —0,53876 0,1203 
38 Н› 29,07 —0,0836 0,1064 0,5752 


Уравнения (1), (3) и (4) описывают изменения концентраций реагентов в процессе, а в системе (2) 
используются мольные доли для описания изменения мольного расхода смеси. В (5) рассматривается 
зависимость изменения температуры реакции от термодинамических характеристик веществ. В (6), (7) 
приведена зависимость для расчета энтальпии образования 1-го компонента и удельной теплоемкости. 

В данной модели прямая задача (решение системы дифференциальных уравнений) является как жесткой, так 
и нежесткой. Метод Гира, используемый для решения модели каталитического риформинга, основан на 
дифференцировании назад [13-15]. Для сравнения результатов и выявления оптимального метода 
использовались явный одношаговый метод численного интегрирования Рунге-Кутты четвертого порядка 
и явный метод численного интегрирования зс1ру.о4еш() языка программирования Руоп. 
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Концепцию автоматизации формирования математического описания кинетики для многостадийных 
химических реакций и численного метода решений проиллюстрируем с помощью методологии ШЕЕО. 

На рис. | представлена контекстная диаграмма процесса автоматизации формирования математического 
описания кинетики для многостадийных химических реакций и решения численным методом. 


Метод формирования 
математического 
описания 


Метод Рунге-Кутты 
4-го порядка 


Метод 
зсру.одешк) Метод Гира 


Стехиометрическая 
а Автоматизация формирования 
пр математического описания кинетики Численное решение 
к я 
для многостадийных химических прямой задачи 
Экспериментальные реакций и численный метод решения 

данные > 

НИ 5 


Разработанный 
программный 
комплекс 


Рис. 1. Контекстная диаграмма «Автоматизация формирования математического описания кинетики для 
многостадийных химических реакций и численный метод решения» 


Входные параметры — стехиометрическая матрица реакции и экспериментальные данные. В их число 
входят начальные данные по концентрации веществ, температуре, предэкспоненциальным множителям, 
энергиям активации. Результатом решения задачи является численное решение прямой задачи. В качестве 
элементов управления выступают численные методы для сравнения результатов, а также метод формирования 
математического описания (СОНДУ). В качестве механизма выступает разработанный программный комплекс 
автоматизации математического описания кинетики многостадийных реакций и решения прямой задачи. На 
рис. 2 приведена декомпозиция процесса автоматизации формирования математического описания кинетики 
для многостадийных химических реакций и решения численным методом на этапы: 

— формирование математического описания; 

— решение прямой задачи. 


Метод Рунге-Кутты 
Метод формирования Метод Гира 4-го порядка 
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матрица 
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Решение прямой Численное решение 
задачи прямой задачи 


Разработанный 
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Рис. 2. Этапы решения задачи «Автоматизация формирования математического 
описания кинетики для многостадийных химических реакций и численный метод решения» 
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Возможности разработанного программного комплекса представлены описанными ниже особенностями. 

При помощи программы можно создать кинетическую модель многостадийной химической реакции 
и представить её решение методами, описанными ранее, используя заданную стехиометрическую матрицу. 
Если реакция неизотермическая и происходит изменение объема реакционной смеси, к системе 
дифференциальных уравнений добавляются два дополнительных уравнения для описания этих изменений. 

Математическая модель каталитического риформинга бензина представляет собой систему 
дифференциальных уравнений (СДУ) вида (1)-(7), состоящую из 40 дифференциальных уравнений по числу 
групповых компонент, температурной зависимости и уравнения изменения молей. 

Основной функционал программы представляет собой: 

— главное окно приложения с выводом математического описания по стехиометрической матрице, 
введенной пользователем; 

— окно с визуализацией решения: для каждого численного метода (Руфоп-одешть, Рунге-Кутта 4-го порядка, 
метод Гира) выводится время работы, массив времени { и решение системы дифференциальных уравнений 
в каждый момент времени, которое представляет собой графики изменения концентрации веществ, изменения 
температуры реакции или объёма реакционной смеси. 

Результаты исследования. Реакция каталитического риформинга протекает в трех реакторах, между 
которыми располагается нагревательная печь. Начальная температура первого реактора равна 766 К (493 °С). 
Первые 9,6 условных часа реакция протекает в первом реакторе, температура падает. После этого происходит 
нагревание в печи между двух реакторов до 763 К (490 °С) и химический процесс переходит во второй реактор. 
Во втором реакторе процесс протекает в интервале времени от 9,6 до 32,3 условных часа, температура 
понижается. Прежде чем перейти к 3 реактору, смесь нагревается во второй промежуточной печи до 768 К 
(495 °С) и протекает в этом реакторе в интервале времени от 32,3 до 60 условных часа. 

На рис. 3—5 изображены кинетические кривые групп углеводородов, участвующие в реакции и полученные с 
помощью метода Гира. Скачкообразные переходы на графиках свидетельствуют о том, что смесь переходит из 
одного реактора в другой с промежуточным нагревом. 
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Рис. 3. Кинетические кривые каталитического риформинга бензина 
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Рис. 4. Кинетические кривые каталитического риформинга бензина 
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Кривые изменения концентраций групповых компонент приведены в заисимости от времени контакта 
реакционной смеси с катализатором (рис. 3—5). Например, символ пР1 обозначает нормальный парафин с одним 
атомом углерода, а символ АЭ — ароматический углеводород с девятью атомами углерода. На рис. б и 7 
представлены соответственно кривые изменения температуры и объема реакционной смеси реакции 


часы 


Рис. 5. Кинетические кривые каталитического риформинга бензина 


каталитического риформинга бензина. 
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Рис. 6. Кривая изменения температуры 
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Рис. 7. Кривая изменения объема реакционной смеси 


Обсуждение и заключение. В представленной работе разработан программный модуль формирования 
математического описания кинетики для сложных многостадийных реакций. Это позволяет анализировать сложные, 
в том числе промышленные процессы, рассматривать разные возможные механизмы за приемлемое время. 

Разработанный программный модуль численного решения прямой кинетической задачи позволяет 
отслеживать динамическое изменение концентраций компонент и температуры процесса в случае его 
неизотермичности. Помимо представленных графиков изменения концентраций компонент, произведены 
замеры времени работы каждого из трех реализованных методов. По методу зсру.одеш() решение прямой 
задачи длилось 80,5 секунды, по методу Рунге-Кутты 4-го порядка — 111,8 секунд. Наиболее оптимальное 
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время, равное 39,8 секунды, показал метод Гира. Это подтверждает, что метод Гира оптимизирован для работы 
с жесткими системами дифференциальных уравнений. 

Реализованный подход применен для моделирования промышленного процесса каталитического 
риформинга бензина. Автоматизация формирования математического описания кинетики и реализация 
численного решения прямой задачи позволили в приемлемые сроки разработать кинетическую модель 
детализированной схемы химических превращений и получить концентрационные и температурные профили 
в реакторах процесса, соответствующие промышленным данным. 
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